Solitons supported by complex PT symmetric Gaussian potentials 
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The existence and stability of fundamental, dipole, and tripole solitons in Kerr nonlinear me- 
dia with parity-time symmetric Gaussian complex potentials are reported. Fundamental solitons 
are stable not only in deep potentials but also in shallow potentials. Dipole and tripole solitons are 
stable only in deep potentials, and tripole solitons are stable in deeper potentials than for dipole soli- 
tons. The stable regions of solitons increase with increasing potential depth. The power of solitons 
increases with increasing propagation constant or decreasing modulation depth of the potentials. 
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I. INTRODUCTION 



Quite recently much attention has been paid to light 
propagation in parity-time (PT) symmetric optical media 
in theory and experiment [Iref. This interest was moti- 
vated by various areas of physics, including quantum field 
theory and mathematical physics [7H15|. Quantum me- 
chanics requires that the spectrum of every physical ob- 
servable quantity is real, thus it must be Hermitian. How- 
ever, Bender et al pointed out that the non-Hermitian 
Hamiltonian with PT symmetry can exhibit entirely real 
spectrum Q. Many people discussed the definition of 
PT symmetric operator and its properties. A Hamilto- 
nian with a complex PT symmetric potential requires 
that the real part of the potential must be an even func- 
tion of position, whereas the imaginary part should be 
odd. It was suggested that in optics the refractive index 
modulation combined with gain and loss regions can play 
a role in the complex PT symmetric potential [16[ . 

Spatial solitons have been studied since their first theo- 
retical prediction [17| ■ Recently, researchers have focused 
on composite multimode solitons. Many composite mul- 
timode solitons are associated with dipole and tripole 
solitons. In local Kerr-type media, fundamental soli- 
tons are stable, whereas multimode solitons are unsta- 
ble. Otherwise, multimode solitons have been studied in 
nonlocal nonlinear media theoretically and experimen- 
tally |18l42(j | . Many authors have paid much attention to 
multimode solitons in optical lattices too |2l| - [23j . 

In this paper, we find that dipole and tripole solitons 
can exist and be stable in Kerr nonlinear media with PT 
symmetric Gaussian complex potentials. The stabilities 
of fundamental, dipole, and tripole solitons are mainly 
determined by their corresponding linear modes for low 
propagation constants or deep potentials. Fundamental 
solitons are stable not only in deep potentials but also 
in shallow potentials. But dipole and tripole solitons are 
only stable in deep potentials, and tripole solitons are 
stable in potentials deeper than that for dipole solitons. 
The stable ranges of solitons increase with increasing po- 
tential depth. 



II. MODEL 

We consider the (1+ l)-dimensional evolution equation 
of beam propagation along the longitudinal direction z 
in Kerr-nonlinear media with complex PT potentials, 

dU d 2 U 
1 -q^ + -q^2+ T[V{x) + iW(x)]U + <j\U\ 2 U = 0. (1) 

Here U is the complex envelop of slowly varying fields, 
x is the transverse coordinate, and z is the propagation 
distance. V(x) and W(x) are the real and imaginary 
parts of the complex potentials, respectively, and T is 
the modulation depth, a = 1 represents the self-focusing 
propagation, and a = represents the linear situation. 
Complex PT symmetric Gaussian potentials are assumed 
as 



V(x) 



W(x) = W xe 



(2) 



where Wq is the strength of the imaginary part. For com- 
plex PT symmetric Gaussian potentials, all eigenvalues 
are real when the real part of the potentials is stronger 
than the imaginary, i.e. Wq < 1.0. Otherwise the eigen- 
values are mixed for Wq > 1.0 [24j |. 

We search for stationary linear modes and soliton so- 
lutions to Eq. |T]) in the form U — f(x) exp(iAz), where 
A is the propagation constant, and f(x) is the complex 
function satisfying the equation 

Xf = -^+T[V(x)+iW(x)]f + a\f\ 2 f. (3) 

We numerically solve Eq. ^ for different parameters by 
the modified square-operator method [251 ] . To examine 
the stability of solitons in PT Gaussian potentials, we 
search for perturbed solutions to Eq. ([TJ in the form 



U = e^{f{x) + [g(x) - h(x)]e dz + [g(x) + h(x)]*e d '}, 

where g(x) <C f(x) and h(x) <C f(x) are the perturba- 
tions, and "*" means complex conjugation. Substituting 
perturbed U(x,z) into Eq. ([T]), linearizing for g(x) and 



h(x), the eigenvalue equations about g{x) and h(x) can 
be derived 
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i[%£ - \h + TVh - iTWg + 2\f\ 2 h 
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(4) 



The growth rate Re(S) can be obtained numerically by 
the original-operator iteration method [2a]. If i?e(<5) > 0, 
solitons are unstable. Otherwise, they are stable. 



III. FUNDAMENTAL SOLITONS 
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We first investigate fundamental solitons in PT sym- 
metric Gaussian potentials with Wo = 0.1. FiguresQJa)- 
QJc) show fields of fundamental solitons with different 
propagation constants and potential depths, which cor- 
respond to the cases represented by circles in Fig. HJa). 
We can see that all the real parts of fields are even sym- 
metric whereas the imaginary parts are odd symmetric. 
With increasing propagation constant A, the beam width 
narrows and the beam intensity increases. With increas- 
ing potential depth, the beam intensity decreases but the 
beam width changes little. 
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FIG. 1. (Color online) Profiles of fundamental solitons with 
Wo = 0.1 at (a) T = 1, A = 0.4; (b)T = 1, A = 2.6; and 
(c)T = 4, A = 2.6. (d)-(f) are the linear modes corresponding 
to (a)-(c), respectively. Solid blue and dotted red lines rep- 
resent the real and imaginary parts of fields, and imaginary 
parts are multiplied by 10. 



FIG. 2. (a) Power P versus propagation constant A with dif- 
ferent modulated depths T for fundamental solitons. Solid 
lines represent stable range and dashed lines represent un- 
stable range, (b Power P versus modulated depth T with 
different propagation constants A for fundamental solitons. 
(c) Perturbation growth rate versus propagation constant A 
with different T values, (d) Critical propagation constant A c 
versus modulated depth T. 



Figures [Ud)-[ljf) are the field distributions of linear 
modes corresponding to Figs. [Ua)-[TIc), respectively. 
The field distributions of linear modes and fundamental 
solitons are homologous for low propagation constants 
[see Figs. [TJa) andQJd)], or f° r deep potentials [Figs. 
[He) andQJf)], but significantly different for large propa- 
gation constants and shallow potentials [Figs. [TJb) and 
[IJe)]. This phenomenon can be explained qualitatively 
by Eq. (fTJ). The nonlinear waveguide produced by the 
term \U\ 2 U, along with the real part of the PT poten- 
tial [V(a;)], confines the expansion of the beam induced 
by diffraction, and also suppresses the transverse energy 
flow induced by the imaginary part of the PT potential 
[W(a;)]. Stationary linear modes or solitons are obtained 
when all these effects are in balance. When the prop- 
agation constant is small, the intensity of fundamental 
solitons and the term \U\ 2 U are small too [see Fig. HJa)]. 
The influence of nonlinearity in Eq. (TTJ) is weak, so the 
field distributions of fundamental solitons are similar to 
those of corresponding linear modes. The field distribu- 
tions of linear modes and fundamental solitons are dif- 
ferent when the nonlinear term is comparable with the 
term V, i.e., for large propagation constants and shallow 
potentials [Figs. OJb) andQJe)]. 

The power of solitons is defined as P — f_ \f(x) \ 2 dx. 
Figure HJa) shows the power of solitons versus the prop- 
agation constant A with different T values, and Fig. &b) 
shows the power of solitons versus the potential depth 
T with different A values. We can see that the power 
of solitons increases with increasing A or decreasing T. 
Figure HJc) shows the perturbation growth rate versus 
the propagation constant A with different T values. One 
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FIG. 3. (Color online) (a)-(c): Evolution of fundamental soli- 
tons corresponding to Figs. QIa)-[l][c), respectively, (d)-(f): 
Evolution of fundamental linear modes corresponding to Figs. 
QJd)-[TJf), respectively. 



can see that the stable range of fundamental solitons is 
A < A c for a fixed T, where A c is a critical propagation 
constant for soliton stability. Figure [U[d) shows that A c 
is approximately proportional to the modulated depth T. 
As T approaches zero, A c approaches zero too, but the 
value of the growth rate Re(8) decreases entirely [see the 
curve of Re(5) for T = 0.01 in Fig. C^c)]. When T = 0, 
Re(6) = in the whole range, and solitons are always 
stable. This is consistent with fundamental solitons in 
pure Kerr nonlinear media always being stable. 

To confirm the results of the linear stability analysis, 
we simulate the soliton propagation based on Eq. ([T]) 
with the input condition U(x,z = 0) = f(x)[l + er]{x)] 
by the split-step Fourier method, where rj(x) is a random 
function with a value of between and 1. e is a pertur- 
bation constant, which is 10 % in our simulation. Figure 
[3] shows the evolution of beams corresponding to those in 
Fig. Q] which is in agreement with the stability analysis 
in Fig. [5J When fundamental solitons are stable, the 
corresponding linear modes propagate with no distortion 
[Figs. 131(d) and|3Jf)]. This means that the linear modes 
absorb the energy of the perturbation noise and main- 
tain its mode profiles. When fundamental solitons are 
unstable, the corresponding linear modes propagate with 
random distortions [Fig. [3Je)]. This indicates that the 
linear modes and perturbation propagate independently 
without energy exchange. According to Figs. [1] and [3J 
we can see that the stability of fundamental solitons are 
mainly decided by the complex potentials for low propa- 
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FIG. 4. (Color online) Profiles of fundamental solitons with 
W = 0.8 at (a) T = 1, A = 0.4; (b)T = 1, A = 2.6; and 
(c)T = 4, A = 2.6. (d)-(f) Linear modes corresponding to 
(a)-(c), respectively. Solid blue and dotted red lines represent 
real and imaginary parts of fields. 

We also study fundamental solitons in the PT symmet- 
ric Gaussian potentials with Wo = 0.8. Figures 0Ja)-|3|c) 
show the field distributions of fundamental solitons with 
different propagation constants and potential depths for 
Wq = 0.8, whereas Figs. |U[d)-|3If) are the linear modes 
corresponding to them. We can see that the properties 
of fundamental solitons for different Wo are very simi- 
lar, except the imaginary parts of fields for Wq — 0.8 are 
larger than those for Wo =0.1. Figure [5] shows the beam 
evolutions corresponding to those in Fig. @] We can see 
that fundamental solitons can propagate stably although 
Wo is close to the point of PT breaking [Figs. [5ja) and 
Uc)}- 



IV. DIPOLE AND TRIPOLE SOLITONS 

We now investigate dipole solitons in PT symmetric 
Gaussian potentials with Wo = 0.1. Figures EJa)-[6jc) 
show the field distributions of dipole solitons, which cor- 
respond to the cases represented by circles in Fig. [TJa). 
Figures EJd)-[S]Jf ) are the linear modes corresponding to 
Figs. E^a)-[6fc) . We can see that all the real parts of 
the fields are odd symmetrical and the imaginary parts 
are even symmetrical, which is converse to the situation 
for fundamental solitons. It is noteworthy that all of the 
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FIG. 5. (Color online) (a)-(c): Evolution of fundamental soli- 
tons corresponding to Figs. |4ja)-|4jc), respectively, (d)-(f): 
Evolution of fundamental linear modes corresponding to Figs. 
ISd)-HIf), respectively. 



field distributions of linear modes and dipole solitons are 
similar in Fig. [6) Due to the deep potentials and small 
propagation constants, the field distributions of solitons 
are decided mainly by PT symmetric Gaussian complex 
potentials. 

The changes in the power versus A and T for dipole soli- 
tons are shown in Figs. [7£a) and[7fb), respectively. The 
power of solitons increases with increasing A or decreas- 
ing T, which is similar to the situation for fundamental 
solitons. Figure [TJ^c) shows the perturbation growth rate 
versus the propagation constant A for different T. Figure 
[7]Jd) shows the critical propagation constant A c versus 
the modulated depth T. We can see that dipole solitons 
exist stably only in deep potentials , i.e. T > 3, and the 
stable range increases with increasing modulation depth 
T. 

Dipole solitons are always unstable for T = 0, which 
corresponds to pure Kerr nonlinear propagation. A 
dipole soliton can be considered two solitons with a it 
phase difference, and a repulsive force exists between 
them. However, we find that the dipole solitons are stable 
in deep PT symmetric Gaussian potentials. The reason 
is that the inherent repulsive interaction between solitons 
can be effectively overcome by the real parts of the com- 
plex PT symmetric potentials. This is the reason that 
dipole solitons are stable only in deep potentials. 

Figure [S] shows the beam evolutions corresponding to 
those in Fig. [BJ which is agreement with the stability 
analysis in Fig. [7J The relations of the propagation be- 
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FIG. 6. (color online) Profiles of dipole solitons with Wo =0.1 
at (a) T = 4, A = 0.3; (b)T = 4, A = 0.8; and (c)T = 
5, A = 0.8. (d)-(f) Linear modes corresponding to (a)-(c), 
respectively. Solid blue and dotted red lines represent real and 
imaginary parts of fields, and imaginary parts are multiplied 
by 2. 
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FIG. 7. (a) Power P versus propagation constant A with 
different modulated depths T for dipole solitons; Solid lines 
represent stable range and dashed lines represent unstable 
range, (b) Power P versus modulated depth T with different 
propagation constants A for dipole solitons. (c) Perturbation 
growth rate versus propagation constant A with different T 
values, (d) Critical propagation constant A c versus modulated 
depth T. 
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FIG. 8. (Color online) (a)-(c): Evolution of dipole solitons 
corresponding to Figs. EJa)-[iHc), respectively, (d)- (f): Evo- 
lution of dipole linear modes corresponding to Figs. EJd)-[6jf), 
respectively. 



tween dipole solitons and the corresponding linear modes 
are similar to those between fundamental solitons and 
their linear modes. 

Finally, we study tripole solitons in PT-invariant Gaus- 
sian potentials with Wo = 0.1. Figure [9] shows the field 
distributions of tripole solitons and their corresponding 
linear modes, which correspond to the cases represented 
by circles in Fig. [TUTa). We can see that all the real 
parts of the fields are even symmetrical and the imagi- 
nary parts are odd symmetrical, which is similar to the 
fundamental solitons. Similar to dipole solitons, tripole 
solitons exist stably in deeper potentials , i.e. T > 8 [Fig. 
[TUT d)]. and all of the field distributions of linear modes 
and dipole solitons are similar. 

The power of solitons increases with increasing A and 
decreasing T, which is similar to fundamental and dipole 
solitons, as shown in Figs. ITUT a) andQjJJb). Figure HUT c) 
shows the perturbation growth rate versus propagation 
constant and Fig. [TOM) shows the critical propagation 
constant A c versus modulated depth. We can see that 
that tripole solitons are stable when T > 8, which is 
larger than the value for dipole solitons. The reason is 
that a tripole soliton can be considered two pairs of out- 
of-phase solitons, and the repulsive force between them is 
stronger than that for dipole solitons. Therefore, it needs 
larger modulation depth T to support tripole solitons 
than that for dipole solitons. 

Figure [TT] shows the beam evolutions corresponding to 
those in Fig. |H1 which are in agreement with the stability 
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FIG. 9. (Color online) Profiles of tripole solitons with Wo = 
0.1 at (a) T = 10, A = 0.3; (b)T = 10, A = 0.6; and (c)T = 
12, A = 0.6. (d)-(f) Linear modes corresponding to (a)-(c), 
respectively. Solid blue and dotted red lines represent real 
and imaginary parts of fields. 
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FIG. 10. (a) Power P versus propagation constant A with 
different modulated depths T for tripole solitons. Solid lines 
represent stable range and dashed lines represent unstable 
range, (b) Power P versus modulated depth T with different 
propagation constants A for tripole solitons. (c) Perturbation 
growth rate versus propagation constant A with different T 
values, (d) Critical propagation constant A c versus modulated 
depth T. 
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FIG. 11. (Color online) (a)-(c): Evolution of tripole solitons 
corresponding to Figs. IHa)-[5Jc), respectively, (d)- (f): Evo- 
lution of tripole linear modes corresponding to Figs. |9jd) 
-IHJf), respectively. 



analysis in Fig. [TU] 



the energy flows from the imaginary part of PT poten- 
tials. Thus fundamental solitons are stable not only in 
deep potentials but also in shallow potentials. For multi- 
pole solitons, a repulsive force exists between their peaks, 
which needs a large modulation depth for counterbalance. 
Therefore multipole solitons exist for a large modulation 
depth and relatively small propagation constant. Our 
results may be extended to other PT symmetric optical 
system, in which multipole solitons can exist. 

Our model [Eqs. (p} and ^] is given in dimension- 
less form, where x and z are scaled to the potential 
width a and the diffraction length L = 2konoa 2 , respec- 
tively. Here kg is the wave number in vacuum and no 
is the refraction index, a is defined as the half width 
at the 1/e maximum of the real part of PT potentials 
[see Eq. © ], so the full width at half- maximum for 
V(x) is about 1.665a and the extreme of W{x) is located 
at ±0.707a. The modulation depth T is scaled to the 
parameter 1 / '{2k 2 iiQa 2 ) . Thus, for a typical waveguide 
a = 10/irn with the substrate refraction index 77o = 3, 
the wavelength Ao = 1.0/jm, and the diffraction length 
L = 3.77mm. Then T = 1 means that the maximum 
variation of refractive index is 4.22 x 10 -5 , and W = 0.1 
means that the maximum gain/loss coefficient is about 
0.53cm . For these physical parameters, it is feasible 
to realize multipole solitons in synthetic PT symmetric 
systems. We hope that the various types of solitons may 
provide alternative methods in potential applications of 
synthetic PT symmetric systems. 



CONCLUSION 
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